Mott transition in anharmonic confinement 
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Two effects are identified that affect the visibility of the Mott transition in an atomic gas in an 
optical lattice confined in a power-law potential. The transition can be made more pronounced 
by increasing the power law, but at the same time, experimental uncertainty in the number of 
particles will induce corresponding fluctuations in the measured condensate fraction. Calculations 
in two dimensions indicate that a potential slightly more flat-bottomed than a quadratic one is to 
be preferred for a wide range of particle number fluctuation size. 



In a spectacular experiment a few years ago, the tran- 
sition from a superfluid to a Mott insulator was realized 
in a gas of cold atoms in an optical lattice by Greiner 
et al. [1]. The experiment marked the birth of a new 
subfield of physics, where theoretical statistical- physics 
models, until then considered to be mostly of academic 
interest or crude approximations to real materials, can 
now be realized in experiment. With the help of atoms 
in optical lattices, it is hoped that fundamental questions 
regarding phase transitions can be addressed, as well as 
applications to e. g. quantum information. Nevertheless, 
a few practical issues remain to be sorted out. One such 
issue concerns the role of the confining potential. 

A gas of spinless bosonic atoms in an optical lattice is 
known to be well described by the Hubbard model [2j, 
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where t is the tunneling matrix element, ^ is the chem- 
ical potential, V(r) is the spatially dependent trapping 
potential, and r is a dimensionless site index. The op- 
erators a r ,al destroy and create a particle at the site 
r, respectively, and obey Bose commutation relations. 
The units are here chosen so that the on-site interaction 
strength, i. e., the prefactor of the first term in ([T]), is 
unity. The sum subscripted < rr' > runs over all pairs of 
nearest-neighbor sites. In the absence of an external po- 
tential V(r), this Hamiltonian exhibits a quantum phase 
transition at zero temperature, separating two ground 
states [|| : When the tunneling matrix element t is strong 
enough, there is phase coherence over the entire sample 
which puts the system in its superfluid state, and there 
exists an accompanying Bose-Einstein condensed fraction 
of the gas that can be measured in time-of-flight experi- 
ments Q. For weaker tunneling, phase coherence is lost, 
the number of atoms per site is locked to an integer, and 
number fluctuations are suppressed. This is the Mott 
insulating state. 

In the experiment reported in Ref. [l|, just as in nearly 
all experiments on optical lattices, the atoms were con- 
tained in a quadratic, i. e., harmonic-oscillator potential, 
in order for them not to escape. As a result, the tran- 
sition between superfluid and Mott insulator is not si- 
multaneous over the whole sample, but takes place via 
an intermediate state where part of the atomic cloud is 



superfluid and part is Mott insulating. Such a state is 
often called a "Mott plateau" state because of the char- 
acteristic density profile, in which the density is fixed to 
an integer in confined regions. This is well known and it 
is also well understood [1,0]. In effect, it is the local con- 
fining potential that contributes a (negative) addition to 
the chemical potential, resulting in a spatially dependent 
critical point for the phase transition. 

In order to see a clear Mott transition one would need 
to get rid of the effects of the confining potential. It has 
therefore been proposed that the harmonic-oscillator po- 
tential could be replaced by a more flat-bottomed variety, 
e. g., a quartic or sixth-power potential, more similar to 
a square well 0]. In such a potential, the local chemical 
potential varies slowly in the center of the sample, where 
most atoms are residing, and the gradual character of the 
Mott transition would be less pronounced. 

However, and this is the theme of the present paper, 
a pure Mott insulating state requires commensurate fill- 
ing, i. e., the number of particles has to match the num- 
ber of wells. In a quadratic potential this requirement is 
relaxed because excess atoms are absorbed into a super- 
fluid region, where the filling is noninteger, at the surface 
of the sample. Using a flat-bottomed potential may put 
more severe constraints on the number of particles, which 
presents a problem, since in an actual experiment one 
does not have very precise control over this number. It is 
the purpose of the present study to explore the balance 
between on the one hand, the gradual character of the 
Mott transition, and on the other hand, the uncertainty 
in the particle number. 

We consider two-dimensional systems in this study, 
since a power-law potential is more easily created in the 
plane than in three dimensions. The confinement to two 
dimensions can be realized by applying a strong optical- 
lattice potential in the third direction. The external po- 
tential is taken to be 



V(r) = (\r\/r o y 
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where |r| is the distance from the center in two dimen- 
sions measured in unit cells, tq is a constant length, and 
p is the power law. In experiments with trapped atoms, 
quadratic potentials are prevailing for natural reasons 
(the first term in the expansion around a minimum of a 
smooth function is in general quadratic). A fourth-order 
potential has been created by superimposing a Gaussian 
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FIG. 1: Cross-section of the density profiles in a two- 
dimensional power-law potential with power law p — 2 (a), 
and p = 8 (b), respectively. The number of atoms is N = 300 
(full lines), N = 400 (dashed), and N = 500 (dotted). 



optical potential on a quadratic magnetic one in order 
to observe fast-rotating vortex configurations [8] . Other 
than that, realizations of high-order power-law poten- 
tials for atoms have been scarce. In two dimensions, 
an optical power-law potential could be realized using 
Gauss-Laguerre optical beams, Fourier optics, or masks; 
a three-dimensional trap would have to be created by 
combining several such potentials. Figure [T] shows a few 
examples of density profiles for a two-dimensional system 
in quadratic (p = 2) and p = 8 power-law potentials. 

The most common way of detecting the state of an 
atomic gas in an optical lattice is by time of flight, that 
is, to release the sample from the trapping potential and 
image it after some time of expansion. From the resulting 
density profile one may deduce the fraction of coherent 
atoms, i. e., the condensate fraction Recently ex- 
perimental techniques have been developed for directly 
detecting Mott plateaus in situ [HEl; nevertheless, that 
type of method addresses the occupation number and 
not the coherence, so the time-of-flight method remains 
the most straightforward. A harmonically trapped gas 
will in such an experiment exhibit a gradually increasing 
condensate fraction as the tunneling is increased. 

In order to observe a sharper Mott transition, it was 
suggested in Ref. Q that the potential can be made more 
similar to a square well. For a power-law potential of the 
form |J3J), the limit p — > oo corresponds to a square well 
of radius tq. In a flat-bottomed potential, the typical sit- 
uation will be that nearly all atoms will be in the same 
state (Mott insulating or super-fluid) , instead of having 
appreciable fractions of atoms in both states for a broad 
range of parameter values. However, now there appears 
a new problem: a flat-bottomed potential is much more 
sensitive to the exact number of particles. In a homo- 
geneous Hubbard model, the Mott insulating state can 
only occur if the filling is commensurate, that is, if the 
number of particles is exactly an integer multiple of the 
number of wells: Adding one particle to a Mott insu- 
lator results in a superfluid. In a power-law potential, 
there is no such restriction on the number of particles, 
but the region in phase space that can accommodate a 
partially Mott insulating system will be smaller for a flat- 



bottomed potential than for a more rounded one. This 
can be seen Fig. [TJ where the density profiles for a one- 
dimensional system with different particle numbers are 
compared. In the quadratic p = 2 trap, an increase of 
the number of particles does not have a dramatic effect 
on the profile or the number of Mott insulating atoms. In 
contrast, for a p = 8 potential illustrated in Fig.Q](b), it 
is seen how an increase in particle number from N = 300 
to 400 takes most atoms out of the Mott phase. 

In order to quantify this, we use the Gutzwiller mean- 
field method with a local-density approximation. The 
fraction of Bose-Einstein condensed atoms is defined as 

Nc = Er\(ar)\ 2 (3) 

N E>W) ' 

This fraction can be inferred from experimental measure- 
ments [4] and is a direct quantitative measure of the state 
of the system. Inverting the definition of a local chemical 
potential, 

M (r)= M -y(r)=M-(^) P , (4) 

one can apply a local-density approximation to transform 
the spatial sum to a sum over \J = /x(r), which in two 
dimensions reads 

Nc = J^%2nr^')n c ^') 
N I^%2^)nW 

where nc(fJ>) — |(a)| 2 and n(/i) = (<ra) are the conden- 
sate density and total density respectively, for a homoge- 
neous sample at chemical potential /i. The lower limit [j,q 
is the chemical potential at which N vanishes; it obeys 
[1q < with the equality holding at vanishing tunnel- 
ing matrix element t. These functions have been com- 
puted and tabulated, whereafter the integrals have been 
summed for a range of particle number N and tunneling 
matrix element t. The result is displayed in Fig. [2] It is 
seen that for more flat-bottomed potentials, the transi- 
tion between the Mott and superfluid states, i. e., from 
zero to non-negligible condensate fraction, is steeper. At 
the same time, the Mott state occupies a much smaller 
region in phase space for the flat-bottomed potentials. It 
is this tradeoff that has to be considered in experiment: 
In order to see a sharp Mott transition, one needs to have 
control over the number of particles. 

Next, we simulate an experiment where the condensate 
fraction is measured for a series of increasing tunneling 
matrix element t (i. e., decreasing optical lattice irradi- 
ance), where N has been allowed to fluctuate. Experi- 
mentally, controlling the number of particles is very dif- 
ficult but the number can be measured afterwards 
to within certain bounds, so on performing a series of 
runs, one may apply post-selection, i. e., only the data 
from those runs that correspond to the desired number 
is included in the analysis. While the relative number 
between different shots is quite easy to measure, infer- 
ring the absolute number is more problematic and may 
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FIG. 2: Phase diagrams for a Bose-Hubbard model in power- 
law trapping potentials with power (a) p = 2, (b) 8, (c) 20, 
and (d) 100. The phase space is spanned by the tunneling 
matrix element t and the number of particles divided by the 
trap area, N/irrg. 



produce a systematic error In the simulation, the 

average particle number (N) is chosen to lie in the most 
pronounced Mott insulating region seen in Fig. [3] Thus, 
for p = 2 we choose (N) = 0.77t?-q, and so on; for p = 20 
we choose (N) = tttq, corresponding to a filling of one 
particle per site. The fluctuations around the mean par- 
ticle number are assumed to be uniform, as they would 
if the number was controlled by post-selection. Figure 
[3] shows the result for a fluctuation of 10 percent. As 
expected, for the quadratic p — 2 potential, the uncer- 
tainty in particle number has little effect, but the transi- 
tion is gradual. In contrast, for the most flat-bottomed 
potential, p = 20, the transition would have been much 
sharper, but the uncertainty in the number of particles 
ruins the observability of this sharp transition. However, 
in the intermediate range of powers, both effects may 
apparently be working in favor: the transition is reason- 
ably sharp as a function of t, yet the Mott insulating 
state covers a broad enough range of particle number N 
that the experimental uncertainty can be accommodated. 
Upon optical inspection, it seems that in this example, 
all power laws p between 2 and 20 yield about equally 
sensible results. 

As an example calculation in order to quantify these 
considerations, we try to locate the transition by fitting 
the simulated data Nc/N to a function of the form 

where t c is a trial transition point and a and (3 are fitting 
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FIG. 3: Outcome of a simulated experiment measuring the 
fraction of condensed atoms Nc/N in a trap, traversing the 
Mott transition as the tunneling t is increased while the num- 
ber of atoms is allowed to fluctuate within 10% according to a 
uniform random distribution. The different curves correspond 
to different power laws p of the two-dimensional trapping po- 
tential. 



parameters. In Fig. [5J the squared residual error of the 
fit. 



A 2 = (\N C /N~f(t)\ 2 ), 



(7) 



is displayed for different power laws p and uncertainties 
in particle number A A. Clearly, the error in the fit is not 
dramatically different for different power laws, and the fi- 
nite number of points included in the statistical sampling 
introduces a scatter in the data, but for very small uncer- 
tainty, the minimum is of course obtained for the most 
flat-bottomed trap, p — 20. For the largest uncertainty 
of 20%, the best fit is obtained for the quartic trap, p = 4, 
where the scatter in the data is smaller. Further numer- 
ical experimentation, not shown here, indicates that the 
p = 2 trap is best when the uncertainty exceeds about 
30%. For the middle-ground cases, with uncertainties of 
5% and 10%, it is an intermediate power law, p = 8, that 
produces the best fit. So indeed, for a modest uncer- 
tainty in particle number there seems to exist a range of 
power law where neither of the two competing processes 
is very strong, and the clearest Mott transition occurs for 
an intermediate power-law potential. 

Summing up, in power-law potentials one faces a trade- 
off between two effects that affect the detectability of a 
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Mott insulator. On the one hand, a flat-bottomed po- 
tential (i.e., a high power law) is preferred in order to 
make the Mott transition sharper. On the other hand, 
such a potential will increase the precision with which the 
number of particles has to be controlled. It was seen in 
an example simulation that for an uncertainty in number 
between ten and twenty percent, a power law between 4 
and 8 yields the most accurate determination of the Mott 
transition; if the uncertainty can be made smaller, then 
a very flat-bottomed potential with a power law above 
10 yields the best result. 
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